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Identifying the time reversal symmetry of spins as a symplectic symmetry, we develop a large N 
approximation for quantum magnetism that embraces both antiferromagnetism and ferromagnetism. 
In SU(N), N > 2, not all spins invert under time reversal, so we have introduced a new large N 
treatment which builds interactions exclusively out of the symplectic subgroup[5P(A'')] of time 
reversing spins, a more stringent condition than the symplectic symmetry of previous SP(N) large 
TV treatments. As a result, we obtain a mean field theory that incorporates the energy cost of 
frustrated bonds. When applied to the frustrated square lattice, the ferromagnetic bonds restore 
the frustration dependence of the critical spin in the Neel phase, and recover the correct frustration 
dependence of the finite temperature Ising transition. 



I. INTRODUCTION 

The search for simple, controlled approximations 
which capture the collective behavior of matter is a key 
goal of condensed matter. In quantum magnetism, this 
search is hindered by the lack of a small parameter; after 
more than a decade, theorists and experimentalists are 
still searching for a physically realizable quantum spin 
liquid 1 , and the ground state behavior of highly frus- 
trated magnets, like the kagome 2 *^, pyrochlore^ and 
hyperkagome^^ lattices is still unclear. One approxima- 
tion that has proven successful is the "large N" expan- 
sion, which generalizes the model of interest to a family of 
models where the number of internal degrees of freedom 
is indexed by an integer N. As N goes to infinity, central 
limit effects permit the underlying collective behavior of 
the model to be solved exactly, and finite TV properties 
may be obtained from a power series expansion in l/N 
about this solution. 

The basic equation of quantum magnetism is the 
Heisenberg Hamiltonian, 

H = ^J i:j S t ■ Sj (1) 

ij 

where the spin on each site, Si lives in the group SU(2). 
The exchange coupling J can be either positive or nega- 
tive, for simple lattices these lead to antiferromagnetic or 
ferromagnetic ground states, respectively. Both ground 
states break both spin rotational and time reversal sym- 
metries, but the antifcrromagnet is invariant under the 
combination of time reversal and translation by one lat- 
tice site. More complicated lattices can lead to spins 
which are not collinera, so called spiral magnets, or pos- 
sibly to a state in which the spins are not ordered at all, 
a spin liquid. 

As we extend the theory of interacting SU(2) spins 
into a family of related theories, we will lose some of the 
physics unique to SU (2) spins, so how do we guaran- 
tee that our resulting theories still capture the defining 
characteristics of magnetism? What defines magnetism? 
What defines a spin? As always, the first, best place 
to look is at the symmetries. An ST/ (2) spin Hamilto- 
nian has two symmetries - time reversal invariance and 



invariance with respect to SU(2) rotations. XY and 
Ising spin Hamiltonians also obey time reversal and ro- 
tational invariance, but under U(l) or Z2 rotations. The 
spins themselves define a unique direction on a mani- 
fold, CP 1 for SU(2), and invert under time reversal, 
S — > — S. The ground state can break the rotational 
and time reversal symmetries, traditionally simultane- 
ously, as in a ferromagnet, but more recently hypothe- 
sized states can break either rotational symmetry, but 
not time inversion; e.g. a spin-nematic defines a unique 
direction, but does not have magnetic long range order—, 
or chiral spin states which break time reversal, but not 
rotational symmetry 10 . Certainly there are nonmagnetic 
states, such as liquid crystal displays, which also break 
rotational symmetry, so the rotational properties of spin 
are not enough to define magnetism. We propose that 
both the rotational and time reversal properties of spins 
are defining symmetries of magnetism, and that a large 
N theory with broad applicability must maintain both of 
these properties in the large N limit. 

Large N theories are well known in particle physics 11 
and in heavy fermion thcor y 12 i 13 ' 14 , but they were first 
introduced in quantum magnetism by Berlin and Kac, 
who solved the spherical model of ferromagnetism ex- 
actly in a large TV limrfei^. Simultaneously, Anderson, 
Dyson and Maleev introduced spin wave theory, which 
takes the spin S to be large and expands in l/ g 16 i 17 i 18 . 
The large S limit is a classical limit, where the spins be- 
have like classical vectors, rotating under the group 0(3). 
The long wavelength fluctuations of quantum spins were 
studied semiclassically in the nonlinear sigma model 19 , 
where quantum renormalizations to the classical param- 
eters were calculated in the large N limit by extending 
the order parameter manifold of SU(2) spins, CP 1 to 
that of SU(N) spins, CP N ~ 1 ^.. 

The large N quantum limit of magnetism was first 
treated by Affleck and Marsto n 21 i 22 . They used a 
fcrmionic spin representation to treat the S = 1/2 
Hubbard and Heisenberg models by extending SU(2) 
to SU(N), preserving the rotational invariance of the 
Hamiltonian under SU (N) . Since then, other extensions 
of SU(2) have been used, including SP(2N) by Ran and 
Wer*22. Fermionic large N theories capture the physics in 
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the extreme quantum limit, S/N <C 1, where the ground 
state is always disordered. These are useful for study- 
ing the possible spin liquid ground states 2 ^, but not for 
determining if a particular model is a spin liquid in the 
first place. For that, one needs a bosonic spin represen- 
tation, where magnetic long range order corresponds to 
the condensation of the bosons. Arovas and Auerbach 
introduced the bosonic SU(N) theory- 5 , which can treat 
arbitrary ratios of S/N, and both magnetically ordered 
and disordered states. This theory was quite successful at 
describing ferromagnets and bipartite antiferromagnets, 
but is unable to treat frustrated magnets. To resolve this, 
Sachdev and Read extended the theory to arbitrary anti- 
ferromagnetic bonds by limiting the rotational invariance 
to the group SP(JV)2&. However, neither of these theories 
preserve the time inversion properties of spins, because 
for N > 2, not all SU(N) spins invert under time re- 
versal, and although Sachdev and Read's Hamiltonian is 
invariant under SP(N) rotations, it still contains SU(N) 
spins with the wrong parity under time reversal. Re- 
cently we introduced a new large N limit, which identifies 
time reversing spins with the generators of SP(N), and 
then builds interactions exclusively from these symplectic 
spins^i. This condition is more stringent than Sachdev 
and Read's, and leads to a unique large N limit which 
we call "symplectic- A" . In collaboration with Dzero, we 
introduced symplectic- N using a fermionic spin repre- 
sentation to treat Kondo physics and superconductivity 
in the two channel Kondo mode&. Here, we develop 
the bosonic symplectic- N approach for the Heisenberg 
model, which enables us to treat ferromagnetism and an- 
tiferromagnetism on equal footing. 

The structure of this paper is as follows. In section 
HI! we show that the time reversal of spins is a sym- 
plectic property, and extend time reversal to large N 
where the SU (N) generators separate into two classes - 
those that reverse under time reversal, and those that do 
nothing. We examine different decouplings of the large 
N Heisenberg hamiltonian and show that excluding the 
non-time reversing spins from the interaction Hamilto- 
nian captures both ferromagnetic and antiferromagnetic 
correlations. In section HITl we derive the mean field equa- 
tions for a generic Heisenberg magnet in the symplectic- 
N limit, while in section [TV] we apply these ideas to the 
two dimensional J\ — J2 model, finding both the zero tem- 
perature and finite temperature phase diagrams. Finally, 
in section fVl we draw conclusions about the application 
of symplectic- N to other models. 



II. TIME REVERSAL AND SYMPLECTIC 
SYMMETRY 

Time reversal is defined by its action on an electron 
wavefunction ?/v (x, t) : 

0Vv(x,t) = aVC(x,-i). (2) 



More generally it is a matrix operator, 9 = eK, where K 
is the complex conjugation operator, Kip — ip* K and e is 
the antisymmetric matrix 102. A consistent definition of 
time reversal requires that 9 commute with the unitary 
rotation operators U, the members of the group SU(2), 

U9U Jf = 9. (3) 

Using the definition of 9 = iK, and noting that K con- 
verts W to U T , we find 

UeU T = e (4) 

This is a symplectic condition on the matrices U because 
it requires the invariance of an antisymmetric matrix e 
under orthogonal transformations. If U is taken to be 
the matrix for an infinitesimal rotation, U = 1 + a ■ S, 
the symplectic condition requires that 

S^dSe- 1 =eS T e = -S. (5) 

So the symplectic condition is equivalent to the time re- 
versal of all ST/ (2) spins. 




FIG. 1: The time reversal operator, 8 and rotation operators 
U acting on a spin must commute, a. Depicts S — » —RS. 

b. S R( — S), where R is the rotation performed by U. 
These two are equivalent for all U in SU(2), and SP(N), but 
not in SU(N). 



A. Time reversal in large TV 

In quantum magnetism, it is convenient to use the 
Schwinger boson representation of spins 25 , 

% = \ b lA h i,f> ( 6 ) 

where bj — (fej +1 , fej _ 1 ) is a two component spinor de- 
fined on each site and S a is one of the SU{2) generators, 
eg - the Pauli matrices. When this treatment is general- 
ized to large N, the number of spin components increases 
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from 2 to an even number TV = 2k. Dropping the site 
index, we have 7? = \^ a T^bj^, where 



(7) 



and T a are the generators of the group SU(N). 

Time reversal is a defining property of magnetism, so 
when we extend the number of spin components, we want 
to maintain this essential discrete symmetry. However, 
SU(N) generators divide into two classes under time re- 
versal(Fig[3D, 



-T a oe{l,2,. 
+T a a e {D N - 



-,D N } 
1,...,7V 2 



1} 



(8) 



where Dm — ~N(N + 1). The first class can be identified 
as the generators of the symplectic subgroup, SP(N), 
whose elements reverse under time reversal, just like the 
SU(2) spins. Dm is the number of N dimensional sym- 
plectic generators. To avoid confusion, we will label these 
symplectic spins by S a . The second class does not in- 
vert under time reversal, and does not form a closed 
subalgebra of SU(N). By analogy with the fermionic 
representation^, we know that these "antisymplectic" 
generators change sign under particle-hole transforma- 
tions, or charge conjugation, while the symplectic spins 
are neutral. These generators then behave like electric 
dipoles, not like magnetic moments, so we label them by 
V a . For N = 2, SU{2) SP{2), and there are no anti- 
symplectic generators. However, for any N > 2, the two 
groups are no longer isomorphic. For example, SU(i) 
consists of ten symplectic generators 



S a e 



il 



-il 



±a 



(9) 



(corresponding to the four Dirac matrices 7 M and their 
six commutators | [7^,71/]), an d five antisymplectic gen- 
erators 



V a e 



1 



1 a 



1 



(10) 



(corresponding to the 75 matrix, and its product with 
the four Dirac matrices 1757,!.) 

Here, the choice of SP(N) is motivated by the de- 
sire to maintain the time reversal symmetry of spin in 
the large N limit, but Sachdev and Read originally de- 
veloped SP(N) because, unlike SU(N) it contains well 
defined particle-particle singlets^. SU(N) expansions 
are extremely useful to particle physicists because there 
are two well defined color singlets - mesons and baryons. 
Mesons are particle-antiparticle pairs, or in condensed 
matter, particle-hole pairs, while baryons are products of 
N particles, forming the three quark baryons for SU(3), 
which have no condensed matter analog except for N = 2, 
where these are particle-particle pairs - e.g. valence 
bonds 28 or Cooper pairs. In the large N limit, the con- 
densed matter version of SU(N) has only particle-hole 



pairs. However, the group SP(N) does have well defined 
particle-particle singlets, which are the pairing of a par- 
ticle and its time reversed twin, and particle-hole pairs, 
but no baryons The presence of these well defined sin- 
glets is equivalent to the existence of a well defined time 
reversal symmetry of spin. 



B. Decouplings of the Heisenberg Hamiltonian 

In order to treat magnetic interactions, we would like 
to rewrite the Heisenberg Hamiltonian, H = JSfSj with- 
out explicit reference to the spin generators. In SU(N), 
this is done by using the SU (N) completeness relation: 



2 A A 



(11) 



We now derive a similar SP(N) completeness relation. 
Any even dimensional matrix can be split into a symplec- 
tic and antisymplectic part: M = Ms + Ma, where the 
symplectic part satisfies Ms = — eMje T and the anti- 



symplectic part Ma — tM A e . The symplectic part can 
be obtained by projection, Ms = PM, where P is defined 
such that PM A = 0. We recognize that M A -iM^i T = 0, 
and take 



PM 



1 



M - eM T e T ) 



(12) 



This expression can be written out in terms of compo- 
nents, 



P^M" 7 = i [M a ^ 

= hs^-e^]M^, 



so that 



pa/3 = _t S a S /3 

iv 2 v 1 



C 7 C 7)J 



(13) 



(14) 



Since the symplectic matrices form a group, Ms can 
always be expanded in the symplectic generators, S a , 
M s = Y,a m ^ a - With the normalization Tr [S a S b ] = 
2S a b, consistent with the SU (2) Pauli matrices, the coeffi- 
cient m a = iTr [S a M], giving PM = | E a Tr [S a M] S a . 
Expanding both sides in terms of components and can- 
celing M 177 , we find 



pct/3 _ 1 \ * ca ca 



(15) 



Finally, by inserting (fT4|) . we obtain the SP(N) com- 
pleteness relation, 



(16) 



Inserting the Schwinger boson spin representation, the 
symplectic N Heisenberg Hamiltonian becomes the sum 
of two terms 



(17) 
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where 

B ii = \ y E 5b i° b $-° ( 18 ) 

cr 

creates a valence bond, or spin singlet pair, between sites 
i and j , and 

4 = ^E^ as) 

a 

creates a ferromagnetic bond, which imply the coher- 
ent hopping of Schwinger bosons from site to site. In 
the language of valence bonds, a ferromagnetic bond can 
be thought of as resonating one end of a valence bond 
between sites i and j, causing both sites to be simula- 
taneously antiferromagnetically correlated with a third 
site, thus ferromagnetically correlated with one another. 
In this sense, it is a frustrating field. Most generally, 
a ferromagnetic bond on a link with antifcrromagnetic 
J, or vice versa, can be considered frustrating fields, 
however, we will usually be dealing with entirely anti- 
ferromagnetic lattices, where any ferromagnetic bond is 
a frustrated bond. This decoupling is identical to the 
££7(2) mean field theory introduced by Ceccatto et al22, 
now controlled by the large N limit of properly time re- 
versing spins. Next, we compare this representation with 
SU{N) 25 and the previous SP(N) treatment. 

The SU(N) Heisenberg Hamiltonian is a dot product 
between SU(N) spins, T, which can be rewritten using 
the SU(N) completeness relation (|TTJ| to obtain the usual 
sum of ferromagnetic bonds 2 ^ 

h suw = J -§%-% (20) 

= ^ (WV\. • S.-S,). (21) 

where Jij is rescaled by N so that H is extensive in N. 
As one would expect in SU (N) , the symplectic and anti- 
symplectic spins are treated on equal footing, which leads 
to a completely ferromagnetic theory. Bipartite antifer- 
romagnets can also be studied in SU(N) by performing 
a special transformation (not time reversal) on one sub- 
lattice, but SU(N) cannot treat more complicated, e.g.- 
frustrated, antiferromagnets. 

The SP(N) Hamiltonian, as defined by Sachdev and 
Read 26 was originally written in terms of valence bonds, 
Hsp(n) = ~JijBj i Bji, in order to treat frustrated anti- 
ferromagnets. When we rewrite it in terms of magnetic 
and electric dipoles, we find 

H SP{ N) = -^B^Bn (22) 
- || (Si'Sj-Vi-Vj) ■ (23) 

Surprisingly, the SP(N) large N theory weights the 
physical symplectic and unphysical antisymplectic spins 



equally, but with opposite signs. SP(N) was so called be- 
cause the Hamiltonian satisfies symplectic symmetry, not 
because it describes the interactions of symplectic spins. 
In fact, any combination of the two terms B^B and A* A 
has symplectic symmetry, including SU(N). The re- 
quirement that our interactions include only mag- 
netic, symplectic spins is more stringent, and this 
method is what we call symplectic- N, while we will 
continue to refer to Sachdev and Read's formulation as 
SP(N). 



Approach 


R(S,V) 


H(6t,6) 


SU(N) 
SP(N) 
Symplectic-N 


J(S-S + P-P) 
J (S -S - P - P) 
JSS 


JA^A 
-JB^B 
J (-B^B + AU) 



Why is it important to exclude the non-time reversing 
dipoles? Both the symplectic {Si ■ Sj ) and antisymplectic 
(Pi ■ Pj) interactions are invariant under time reversal, 
however, the important difference is not in the Hamilto- 
nian, but in the ground states, and the dynamics. These 
are far more coupled than the Hamiltonian suggests be- 
cause the SU(N) spin T does not act as a vector, and the 
antisymplectic and symplectic directions are not indepen- 
dent directions, so that T is unable to point in a purely 
symplectic direction. The antisymplectic interactions en- 
courage the antisymplectic spins to order - competing 
with the ordering of the physical components. This com- 
petition eliminates the antiferromagnetic [ferromagnetic] 
ground state completely for SU(N)[SP(N)}. And finally, 
even if the ground state is the one of interest, the pres- 
ence of antisymplectic interactions affects the dynamics 
of the symplectic spins, dynamically violating the closure 
of the symplectic subgroup. 

C. Constraints and spin Casimirs 

In the Schwinger boson representation, the total spin 
on site is free to take any value. Therefore, in order 
to treat a real spin problem, the total spin must be re- 
stricted to its physical value, S 2 = S(S + 1). This con- 
traint is implemented by a Lagrange multiplier fixing 
the value of the spin Casimir (Sj) 2 , which depends on 
the group. For a general group with generators r a , the 
Casimir is written, 

= EC^OMC^r^)- (24) 

a 

For symplectic spins, the completeness relation (|16j) is 
used to rewrite the Casimir as 

= \ ( b ]a b jp) ( b j-y b 3l) [Sa V S^ + 6^6^] 

= 4 {'^J'j A A/- + & McPj-P h \-J>jp) 
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b ]a b 3* 



t 

30 



n bj 



where n bj = £ Q bj a bj 
j. The last equality is due to the vanishing of antisym- 
metric combinations of bosons, 
for symplectic N, the Casimir is given by 



is the number of bosons on a site 



on site. Thus, 



1 



Without loss of generality we can choose x to be the spin 
-defined by b T = (1,0, .. .0). Rotating b by any matrix 
(25^hich affects only the lowest N — 2 entries clearly leaves 
x invariant, as does rotating the phase of the upper two 
entries, so H x = SU(N - 2) x [7(1), and 



(26) 



M SU (N) = SU(N)/SU(N - 2)x[/(l) £ CP 



N-l 



(30) 



The full SU(N) spin lives on the manifold CP"- 1 , while 
the symplectic spin jbyS-bj lives on a 2N— 1 dimensional 
manifold, given by 



and is set by fixing the number of bosons on each site. 
If we choose the convention n b j = NS, the constraint 
becomes 



M SP(N) = SP(N)/SP(N-2)xU(l). 



(31) 



S? = \n 2 S{S + 1) 



The Casimir for SU(N) can be obtained similarly, us- 
ing the SU(N) completeness relation (TTT]) instead of ([Tfi)) : 



AT 



(27) 



where we have dropped the j index on n b for clarity. 
Using the consistent convention n b = NS, the SU(N) 
constraint becomes 

T 2 = l(N 2 ~N)S(S + l). 
J 2 

For N — 2, this reduces to 5(5 + 1) and the SU(N) and 
SP(N) Casimirs are identical, as required. For all other 
N , T? will be larger. This means that the antisymplectic 

spins, V 2 = 7) 2 — Sj can never be removed for any N > 
2. In the large N limit, they are forced to have equal 
magnitudes: Vj = Sj. 

At first sight, this requirement is quite strange. Af- 
ter all, there are TV 2 — 1 independent SU (N) generators, 
which we have been treating as a vector, T, why can the 
spin not point in A^ 2 — 1 directions? The answer is that 
not all directions of the SU(N) vector give rise to differ- 
ent spins. The spin itself is given by ■ T ■ bj , and b has 
N components. The constraint removes one more degree 
of freedom. For a general state, & is a bosonic vector, but 
when the spins order, 

(4) = \{b) a S a0 b j0 ) = \{tf) 3a S aP {b) 3 p, (28) 

(6) is an N component complex vector, so the spin can 
only take on 2A^ — 1 different configurations. The spins 
are constrained to a 2 N — 1 dimensional manifold M. . 

To be more mathematically precise, this manifold M. 
is a "homogeneous space" of SU(N): SU(N)/H X , where 
H x is the "stabilizer" of x, the subgroup which leaves an 
SU(N) element x invariant: 



Since V is nonzero, M. spin) is not contained within 
Msp(n)] m fact, the two manifolds have equal dimen- 
sion, although they are not isomorphic. Rather, any 
point on Msp(n) wm correspond to a point on Msu(N)- 
Strictly speaking, this manifold is the order parameter 
manifold for a long range ordered state, however, it paints 
a useful picture of the relationship between SU(N) and 
SP(N) spins. Furthermore, the order parameter man- 
ifold will be essential in describing the ordered state, 
where, for a spiral state which completely breaks the sym- 
metry, the number of Goldstone modes will be 2 N — 1 . 



D. Ground States 

A generic Heisenberg Hamiltonian with symplectic in- 
variance contains both antisymplectic and symplectic in- 
teraction terms, 



H = ^./,,5, -5, • i\,,v,- r, 



in a ratio K/J, which is ±1 for SU(N) and SP{N), re- 
spectively, and zero for symplectic- AT. In general, the 
physical, symplectic spins and and the antisymplectic 
spins may have different interaction strengths and signs. 

In the 5 — > oo classical limit, the system is long range 
ordered, and all the bosons are condensed. The ordered 
state is described by the angle between neighboring spins, 
4>ij = <pi — <j)j which is for a ferromagnet and 7r for an 
antifcrromagnet. If we fix (b)i = \A/V5(1, 0, . . .) T , we can 
rotate the top two coordinates of (b)j by 

R(<fo)=( cosi L sin S~ I. (33) 



which makes (b)j = V NS (sin ^ , cos ^ , 0, . . . ) , and 
the two bond expectation values will be 



H x = {g&SU{N)\g-x = x}. 



(29) 



B u = (b)J{b)i=NSsin^ 3 
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Aii = (b)f (b)^NS cos (34) 
Thus the ground state energy for ([32]) is 

E = J2i K - J)ij sin2 ^ + ( K + J)a c o s2 ^ (35) 

ij 

The three special cases of interest are 

Esymp-N = ^ NS 2 Jij COS (j>ij 
ij 

E SP (N) =^-iV5 2 J ii sin 2 ^ 

ij 

E SU(N) = NS2J V C0S2 ^f- ( 36 ) 

ij 

We see that for antiferromagnetic bonds in SU(N), the 
ground state energy is zero, identical to that of the para- 
magnet with (S) = 0, and similarly for the ferromagnetic 
bonds in SP(N). Only symplectic- has well defined 
ground states for both signs of J. 







b. 







FIG. 2: (Color online) A toy picture of SU(N) spins, where 
the symplectic, time reversing components are represented in 
blue along the y axis, and the antisymplectic, non time revers- 
ing components are in red, along the x axis. The purple spin 
shows the full SU(N) spin obtained by adding its symplectic 
and antisymplectic components, (a) depicts a ferromagnetic 
state, (b) depicts an antiferromagnetic state, where we obtain 
the antiferromagnet by time reversing every other spin. While 
the symplectic components are antiparallel, the antisymplec- 
tic components are still aligned, causing the total spins to be 
orthogonal at neighboring sites. 

If we turn to hnitc S, we can construct ferromagnetic 
and antiferromagnetic states explicitly out of the SU (N) 
spins, see Fig [2] The antiferromagnetic state is defined 
by dividing the spins into two ferromagnetic sublattices, 
where sublattice B spins are the time reverse of A. This 
state satisfies the lattice translation plus time reversal 
symmetry of the SU(2) antiferromagnetic ground state. 
The V's are aligned in both ground states, and in the 
large N limit, the magnitudes of S and V are the same. 
In SU(N), K — J, so both interactions are maximally 
satisfied in the ferromagnet - leading in fact to overstabi- 
lization due to excess V bonds, while the antiferromagnet 



consists of orthogonal SU (N) spins, the two terms in the 
Hamiltonian cancel and the ground state energy is zero, 
just as found classically. When K = —J, as in SP(N), 
it is the antiferromagnetic ground state that is oversta- 
bilizcd by ferromagnetic V bonds, and the ferromagnetic 
state has zero energy. These conclusions hold not only for 
the full ground state, but for individual bonds; in frus- 
trated lattices there will be both antiferromagnetic and 
ferromagnetic correlations, even if all J's are positive, 
but SP(N) indicates ferromagnetic correlations only by 
the absence of a bond. The energy cost of ferromagnetic 
correlations is zero in SP(N), but we know that in real 
lattices these frustrated bonds carry a price. By elimi- 
nating the antisymplectic interactions, symplectic- N re- 
moves the extraneous bonds between P's and restores the 
ability of SU(2) to simultaneously treat both ferromag- 
netism and antiferromagnetism. 



E. Spin dynamics 

Even if the ground state is correct, as for the bipar- 
tite antiferromagnet in SP(N), we still need to be con- 
cerned about the spin dynamics. We chose to use the 
group SP(N) not only because its spins all invert un- 
der time reversal, but because the group contains well 
defined particle-particle singlets. The presence of anti- 
symplectic interactions, even if they are only interacting 
with themselves dynamically violates the closure of the 
symplectic subgroup. 

The dynamics of a symplectic spin component at a site 
i are given by 



kj 



K, 



Si,Pk-Vj 



(37) 

We concentrate on the effect of the second term, which 
is generally nonzero when K is nonzero. Inserting the 
Schwinger boson representation, we find 



1 



= kb\-V b -b 3 , b\-S a -b^bl-V b -b k } (38) 



b\-S a -b u 



bl-P" 



b k 



b\-V b -h 



8 l 



where {, } denotes the anticommutator. Expanding out 
the commutator in more detail, 



b]-S a -h,b 



V b -b k 



SapP\r] 



b l a bi/3,bl x b kri 



= 8 ik b\ a [S\V b ] a pb ifi 
= 2i5 lk g ab c t-- (39) 



where g c is the appropriate SU{N) structure factor. 
Since the commutator, [S,V] is odd under time reversal, 
must be an antisymplectic spin. So the evolution of 
Si is affected by the antisymplectic spins, 




v-v 



3 



(40) 
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where x is the cross product defined by g ab c . The full 
dynamics of the symplectic spins are given by 



dSj 
dt 



— [ Si x ^ JijSj 



(41) 



These dynamics are identical in form to classical spin 
wave theory where the spins are torqued by an effective 
magnetic field coming from neighboring spins. The sym- 
plectic and antisymplectic components of % are torqued 
by the effective magnetic fields given by ^ • JijSj and 

"YlijKijVj, respectively. The effective field coming from 
the antisymplectic components is not strictly a magnetic 
field, as it has even time reversal parity, but most im- 
portantly, it rotates V into iS, and vice versa. Ordinary 
SU(2) spin waves will also break spin singlets, but the 
excitations remain in the SU (2) space, while the excita- 
tions for K 7^ will take us out of the SP(N) group, 
ft is clear that to have a theory of interacting SP(N) 
spins, all the antisymplectic interactions must be elimi- 
nated; all other Hamiltonians with symplectic invariancc 
describe anisotropic SU{N) spin interactions. 




SU(N) spins 



H[S] 



p LANE 



FIG. 3: SU(N) spins consist of two components: symplectic 
directions that reverse under time reversal, and antisymplec- 
tic directions that are invariant under time reversal, which 
prevent SU(N) spins from forming two particle singlets. How- 
ever, if the spins are projected into the symplectic plane, these 
components can form two particle singlets, which are well de- 
fined as long as the antisymplectic components are noninter- 
acting. 

So we have seen that the inclusion of antisymplec- 
tic spin interactions have a rather serious effect on the 
physics of the Heisenberg model. When these interac- 
tions are excluded, as in symplectic-iV, the unphysical 
antisymplectic spins can no longer affect the physical 
spins. In a sense, they come along for the ride, since 



they are always there, and they are affected by the sym- 
plectic spins, but have no effect on the physics. Now we 
move on to the application of symplectic- N to general 
lattices. 



III. SOLVING THE SYMPLECTIC- AT 
HEISENBERG MODEL 

Now we return to the symplectic- TV Heisenberg 
model (fl"7|) to discuss how to solve the Hamiltonian in 
the large N limit, for a general lattice specified by Jij. 
As a refresher, the Hamiltonian is 



'"> = E ^ • 4 = E t [-4^ + 4^ 



(42) 



and A] 



hb]„bj a and the sum 

2 l(J J 



where Bj 
over <7 is implied 

The usual prescription for solving these problems is to 
write the partition function as a path integral, 



Z = J Vbc- NS W ]js (S](r) - N 2 S(S + 1; 



where NS[b] is the action 



NS[b] 



(It 



(43) 



(44) 



and the constraint ]\ j 8 (s]{t) - N 2 S(S + 1)) restricts 

the spins to the physical subspace at every site j and time 
r. This constraint can be rewritten using a Lagrange 
multiplier Aj(r), 



(45) 



From now on we drop the explicit r dependence of 
and \i. 

In order to evaluate the path integral, Z must be in 
the form of a Gaussian integral, so the quartic terms in H 
are decoupled using the Hubbard-Stratonovich identity, 



N 
2ttU 



VAc -NAA/J = i 



After inserting this identity, A can be shifted to A - 
eliminating the quartic term j^BB, 



e« BB oc 



DA e - NAA /J+ AB + B &. 



(46) 



(47) 



Now we have exchanged a theory of bosons with four 
particle interactions for a theory of free bosons interact- 
ing with a fluctuating field A. We can integrate out the 
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bosons exactly, but we will need to use the saddle point 
approximation to perform the path integral over A (Fig 
IDJa)), an approximation that becomes exact in the large 
TV limit due to the extensive dependence of the action 
NS on N. First we must treat the other quartic term, 
— jjAA. Naively, we would just change the sign in the 
exponentional in (|47p . which gives 



Aij Aij 



AA cx / x>he +N ' lh/J ~ JlA ~ Ah 



(48) 



However, we must be careful, as the quadratic h term now 
has a positive sign, and the path integral over h appears 
not to converge. To understand this, we step back to a 
simpler case, where A is real and we decouple it with the 
real field a. We begin with e~ Na I J , and can rewrite 

-Na 2 /J = +N(ia) 2 /J ^ +(ia + —A) 2 /J 

= -Na 2 /J + 2iAa+ ^A 2 , (49) 

so that the quartic term — j?A 2 becomes —Na 2 /J+2iAa. 
We now define the mean field value of ia = Hq to be real. 
In fact, let's redefine ia = h = Iiq + iSa, and the identity 
becomes 

e-^ 2 = Jvhe Nh ^ J+2Ah , (50) 
which holds as long as h is integrated along the imag- 





y 


V 




A=x+iy 


1 h=u + iv 




X 








r 



a. 



FIG. 4: Integrating out the fluctuations, (a.) A is integrated 
along real axes x and y, and its saddle point is a minima at 
some Ao. (b.) h is integrated along imaginary axes u and v, 
with a maximum, real saddle point ho- 

inary axis, with the integral maximized at a real ho(see 
FigHJb)). This can be generalized to a complex a = 
u + iv, where u, v are imaginary instead of real. As long 
as we keep in mind that A is integrated along the real 
axis and h along the imaginary axis, we can proceed with 
the above decouplings, 



H[b] = £(«) (6* ah 



A, 




■J, 



(51) 



where y] 
nonzero J;. 



(y) 



is performed only over bonds (ij) with 
The notation can be simplified by defin- 
ing the Nambu spinor, bj = (bj a ,abj- a Y We now have 



the partition function 

Z = J T> [b, A, h, X] e~ 



NS[b,A,h,\] 



(52) 



where the action can be compactly written 

" 1? 



-bi (iuj n T 3 + G;/) bj 



NS[b,A,h,X}= J2 

+ -j-(hijhij - A i3 A lj ) + i\ i N(S + -)d i3 



p— l _ | iXiSij — 2hij 

10 ~ 2A,;,- 



2A y - _ 



(53) 



We have performed a Fourier transform in imaginary 
time, and bi, hij, A^- and Aj are now functions of the Mat- 
subara frequencies iui n , although in practice we make the 
Ansatz that hij, A^- and Xi are all static quantities. The 
factors of i come from rewriting Xibiabiu in terms of the 
Nambu spinors, 6j. 

We can calculate the mean field values of 
, Aij by approximating the path integral 
J V [A, h, X] exp (—NS[b, A, h, A]) by its saddle point 
value, which becomes exact in the large N limit. By 
minimizing the action with respect to hij, Aij and A^, 
we find 



, Hj — 

A,, = 
NS = 



^(btM 

(blh a ) 



(54) 



where (• • • ) denotes the thermal expectation value. How- 
ever, it is simpler to eliminate the bosons altogether by 
integrating them out. 

In order to proceed further, we must make some Ansatz 
about hij, A^ and A^. In principle, hij and Ay can take 
different values on every bond, but for spatially uniform 
states, we choose an Ansatz with the unit cell of the 
lattice, where h and A are defined for each different Jy. 
If Jij — on any bond, so must hij and Ay. We make 
the approximation that iXi(r) = A on every site, taking a 
local constraint and enforcing it only globally. As usual, 
this approximation becomes exact in the large N limit. 

For a square lattice with only nearest neighbor cou- 
plings, this leads to three parameters, which can be fur- 
ther simplified to just A and A, as there are no frustrat- 
ing interactions. Thus we recover the unfrustrated square 
lattice as previously studied in SP(N} 26 . However, for 
frustrated lattices we cannot generally exclude either h 
or A. 
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Sometimes the uniform state will not be sufficient. The 
ground state might break lattice rotational symmetries, 
in which case A^+j and Aj^+g will be different, or trans- 
lational symmetry, requiring Ai_i +X ^ &i+x,i+2£- When 
rotational symmetry is broken, A will remain the same 
on every site, but broken translation symmetry requires 
Xi =/= Xi+x- Since the unit cell is enlarged, there will be 
more than one branch of oj k , which must be summed over. 
However, as long as the state may be specified by a fi- 
nite number of parameters, it may be modeled within 
symplectic-A^. Problems with infinite parameter sets, 
e.g.- spin glasses^, can also be treated within symplectic- 
N , but require more complicated theoretical machinery, 
and will not be treated here. For the rest of this paper, 
we assume translational symmetry, with i\i = A, but this 
treatment can be easily generalized. 

The Fourier transform of the bosonic Hamiltonian, 



G" 1 is 



A - 2h k 2A k 
2A k X-2hi 



(55) 



We can now perform a Bogoliubov transformation 



det(wT3 — G, 



to obtain 



4A L 



(56) 



and integrate out the bosons to obtain the free energy, 
F[h, A, A] = — /3 -1 Tr log Z[b, h,A,X], where the trace is 
over sites and the Matsubara frequencies iu) n , in 

addition to the bosonic degrees of freedom. 



F = A^^log 



2 sinh 



N 



E 7- ( A « A « 



(>:,) 



hijhij) - XNAf s (S + -)(57) 



Let us say we have a set of {hi, h%, ■ ■ •} and 
{Ai, A 2 , . . .}, which have the Fourier transforms 



hk 
A fc 



^ h a J ak 
a 

^ A a S ak 



(58) 



where a labels a bond. The symmetry properties of hij = 
hji and Ay- = — Aji force j a k and S ak to be symmetric 
and antisymmetric in k, respectively. The free energy is 
now 



NAf s 



2 sinh 



J2y a (\&a\ 2 -\h a \ 2 )-X(S+l) (59) 



where z a is the number of bonds of type a per unit cell 
- for a simple square lattice this is just the coordination 



number z = 4. The free energy is now minimized by 
solving the mean field equations dF/dX, dF/dh ai and 
8F/dA a : 



(60) 



_ 1 1 x ^X-2h k 
2z a h a 1 (A - 2h k ) 2j ak ( \ \ 



Ja 



k 



2z a A a _ 1 x - 2A k S ak ( 1 

Ja M s ^ LU k [ nk+ 2 

n k is the Bose function (eP Uh — l) . 



(62) 



iu k = J[X-2h(i 



A. Simple Example 

Now we examine a simple model in detail, the two 
dimensional bipartite square lattice. We know the mean 
field value of h must be zero, however, for pedagogical 
purposes we keep both h and A. 

cosfcy)] 2 — 4A 2 (sinfe E + sinfe^) 2 

(63) 

We wish to minimize the free energy, however, we must 
be careful because h and A are integrated along the imag- 
inary axis. In fact, the free energy should be maximized 
along h and A directions and minimized along A. To ex- 
amine the nature of the extremum, we look at the Hessian 

/ 

i>.\- ,l.\<)h <>\,IA 1 

(64) 



X 



dXdh 
d 2 F 



d 2 F_ 

dh 2 
d 2 F 



dW dXdh OXdA 
pi 2 1? n 2 p O F 



\ dXdA dhdA dA 2 

where A and h are both zero, which is the global mini- 
mum if the temperature is well above where A acquires 
an expectation value. All off diagonal terms vanish at 
this point, 



X = 



/_i csc h 2 f 








V 



_8 _/3csch 2 ^ 

!) 

(65) 

Looking at A and h independently, F is always maxi- 
mized, as expected at the mean field values of h and 
A, while F is minimized along A for small J and maxi- 
mized for large J, indicating a second order transition to 
nonzero A at some intermediate J, dependent on tem- 
perature and spin. 



B. Examining the Ground State 

At zero temperature, we are interested in the ground 
state energy, 

^ = ^5>+Ej(i^-w°)-^4>, 

k a 

(66) 
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which must again be minimized with respect to the pa- 
rameters A, h a , and A a . The order of limits is important; 
to obtain the correct mean field equations or we must 
take the derivatives of the free energy first and then take 
the limit T — * 0. In the mean field equations (|60l - [62]) . 
all temperature dependence is in n k - If there is no long 
range order, limy— >o n k = 0. The Mermin- Wagner theo- 
rem forbids the breaking of a continuous symmetry, like 
SU(2) or SP(N) at any finite temperature in one and two 
dimensions 31 , however at T = 0, the Heisenberg mag- 
net may develop long range order, which corresponds to 
the condensation of the Schwinger bosons^. The bosons 
themselves develop an expectation value, 



+ 2 E n ^/ 2 (A-2^ /2 ) (71) 
Q 

k Q 

Now we have set up all the machinery for solving the 
symplectic-iV Heisenberg model on a general one or two 
dimensional lattice(three dimensional lattices cannot cur- 
rently be treated by Schwinger bosons^). Next we treat 
a simple example which highlights the differences be- 
tween symplectic-./V and previous large N treatments, 
the Ji — J2 model. 



b ia = {b) i + 5b ia , (67) 

where (b)i and (b)i are no longer independent variables. 
Instead (6)j is a complex N component vector, and (b)i = 
. n k will no longer vanish for all k. To see the effects 
of the long range order, we examine the action (|53p again, 
inserting ([67]) 

NS[b, A, h, A] - fdLoJ2 \{b)\ (ium + G^ 1 ) (6), + NS[5b] 

(ij) 

Q 

The linear terms proportional to Sb must vanish, and so 
have been neglected. Q/2 are the zeroes of the Schwinger 
boson spectrum. For ferromagnetism, Q/2 = (0,0), while 
for antiferromagnetism, Q/2 — (ir/2,ir/2). The long 
range order is indicated by the ordering of the spins, 
which are the combination of two Schwinger bosons, so 
the Goldstone modes in classical spin wave theory will 
be given by Q/2 ± Q /2 = and Q, which gives the 
traditional (71", it) ordering vector for antiferromagnetism. 
Now, in addition to the mean field equations, we have the 
condition 

dS/d{b)^ = G Q ^ 2 (6) Q> = 

= ^Q/2< & >Q/2=0. (69) 

So either (6)^ 2 = 0, and we proceed as before, or ^q/ 2 = 
0, which allows us to find the value of {b) qj 2 m addition 
to the original parameters. In fact, = ^13/2(^)^/2' 

so we can simply define n — n Q/2/ LU Q/2 an< ^ ^he mean 
field equations become 

k Q 
1 ^ (2h k - A) 7 afe 



IV. ILLUSTRATION: Ji - J 2 MODEL 

The J\ — J 2 Heisenberg model is one of the simplest 
two dimensional frustrated magnets, 

H = J\ ^ S x ■ Sx+fi + J-2 ^ S x -S x+f2 >, (73) 

where J\ and J 2 describe nearest and next nearest neigh- 
bor interactions, respectively. We consider only antifcr- 
romagnetic J\ and J 2 . 




FIG. 5: (Color online)The Ji — J2 model, (a) depicts an- 
tiferromagnetic order as described by antiferromagnetic va- 
lence bonds(blue) and ferromagnetic bonds(red, dashed), and 
the spin order Q = (h,tt). (b) depicts the collinear order, 
Q — (0,7r), where there are two different antiferromagnetic 
valence bonds(blue and green) and ferromagnetic bonds(red, 
dashed) . 

For Ji ^> J2, the ground state is a Neel antiferromag- 
net, with Q = (it, n) long range order, as long as the spin 
S is greater than a critical spin S c ~ .2. The next nearest 
neighbors are ferromagnetically aligned, so J 2 introduces 
frustration which begins to surpress long range order by 
increasing the critical spin. 
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For J2 3> Ji , the classical ground state consists of two 
interpenetrating but decoupled Neel sublattices. For any 
finite J\, both quantum and thermal fluctuations couple 
the sublattices together through the process of "order 
from disorder," 34 i 35 ' 36 which leads to a long range or- 
dered state with Q = (0, tt) or (it, 0). This transition 
spontaneously breaks the Z4 lattice symmetry down to 
Z2, which, as a Ising symmetry breaking, can survive 
to finite temperatures, despite the loss of the underly- 
ing long range magnetic order 37 . In real materials, this 
transition couples to the lattice and causes a structural 
transition from tetragonal to orthorhomic symmetry 3 ^. 

The phase boundary between the two classical ground 
states is at J\ = 2J2. Conventional spin wave the- 
ory predicts that the ordered moments of both states 
are suppressed to zero even for S — > 00 at this critical 
point, leaving a quantum spin liquid state that exists for 
a small, but finite range of J2/J1 for the physical spin 
S = 1/2 39 . However, at this point, the 1/5 expansion 
fails 40 , and much more theoretical work has been done 
to see if quantum fluctuations stabilize or destabilize the 
spin liquid regio n 41 i 42 ' 43 . The current consensus is that 
the spin liquid ground state is most likely stable between 
■4 < J2/J1 < .6 for S =1/2. 

The J\ — J2 model is an ideal demonstration of the 
importance of ferromagnetic bonds because they enforce 
the frustration price in both Neel and collinear phases. 
This is most obvious on the Neel side, where, without 
ferromagnetic bonds, the state remains unchanged as J2 
increases, until a first order transition to the collinear 
state. The ferromagnetic bonds also enable us to obtain 
the correct temperature dependence of the Ising transi- 
tion temperature in a large N theory. 

A. Valence bond structure 

First we need to describe the relevant states within our 
valence bond picture. We assign a A to each antiferro- 
magnetic bond and an h to each ferromagnetic bond. On 
the Neel side, we have A on all nearest neighbor bonds 
and hd on the frustrating diagonal bonds, as shown in 
Figure [5](a) . This leads to the dispersion relation: 

u% = ^{\-±h d c x c v y-A^{ Sx + Sy y). (74) 

In the collinear state, we must allow the breaking of 
lattice symmetry and consider both h x , h y and A x , A y 
on the nearest neighbor bonds, and hd, A d on diagonal 
bonds. In fact, there are two distinct diagonal bonds cor- 
responding to what would be the two decoupled sublat- 
tices (see Fig(5Jb)). Their magnitude must be the same, 
but the phase between them leads to a U(l) gauge sym- 
metry. If we fix the phase to be tt it is most natural to 
break the lattice symmetry explicitly^ 4 - and choose only 
h x and A y to be nonzero of the nearest neighbor bonds, 
and hd = 0, which gives the dispersion 

oJ c k = y/(X- 2h xCx ) 2 - (AA d c x s y + 2A ySy Y). (75) 



B. T=0 phase diagram 

To examine the frustrating effects of the ferromagnetic 
bonds, we focus on the border between long and short 
range orders at T = 0, as a function of spin, S = n^/N 
and frustration, J%/Ji . The more stable the phase, the 
larger the region of long range order. Long range order is 
lost as the spin decreases below a critical spin S c , which 
is approximately 1/5 for the unfrustrated Neel lattice, 
e.g. both J2/J1 = and J%/J\ = 00. To compare our 
results to the original SP(N)^, we calculate the phase 
boundaries both with h free and with h set to zero. 

Since we are interested in S c , the onset of long range 
order, we know that both n = and <-^q/ 2 = ^- The 

Schwinger boson gap is at Q/2 = (ir/2, ir/2) for the Neel 
phase, and Q/2 = (0, 7r/2) for the collinear state. The 
gap equation, Wq/ 2 = and the mean field equations for 
/i (|7Tj) and A(|72[) can be used to solve for the mean field 
parameters, and then equation (|70|) defines S c , 

1 f \-2h k 
Sc+ 2^] k -^-- (?6) 

Results from these calculations for both symplectic- 
N and SP(N) are shown in Fig [5J For comparison, we 
have also drawn the phase boundaries given by conven- 
tional spin wave theory 3 ^. First let us discuss the re- 
sults far from the critical value of frustration J\ w 2J2. 
The results are most dramatic for the Neel state, where 
SP(N) is oblivious to the frustrating effects of the diag- 
onal bonds, drastically overestimating the critical spin. 
For the collinear state, SP(N) neglects the frustrating 
h x , again overestimating the stability of the long range 
ordered state. On the other hand, symplectic-./V tracks 
conventional spin wave theory for small amounts of frus- 
tration, but they differ in a wide range around the criti- 
cal J2/J1, where conventional spin wave theory is known 
to fail, and a spin liquid ground state is predicted for 
S = 1/2. For symplectic-iV, we calculated the location 
of the first order transition between Neel and collinear 
long range orders by comparing the ground state ener- 
gies of both states(see Appendix Symplectic-iV indi- 
cates a weakly first order transition for S = 1/2, with no 
intervening quantum spin liquid, however, 1/N correc- 
tions, calculated as Gaussian fluctuations from Ceccatto 
et al.'s mean field theory lead to a small region of spin 
liquid for .53 < J2/J1 < ■64 4 i. 



C. Finite Temperatures: the Ising transition 

Now we turn our focus to the J2 ^> Ji side of the phase 
transition and examine the finite temperature Ising tran- 
sition between decoupled sublattices and the collinear 
phase. This phase transition has several possible exper- 
imental realizations, most prominently and recently in 
the iron arsenides 46 ' 47 '— i 49 i 50 . 
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Spin Waves 



eigenvalues of the Hessian of the free energy, 




FIG. 6: (Color online)We compare the critical spin S c = 
(tt) , below which there is no long range order in the ground 
state, calculated within SP(N) (bold red line), symplectic- N 
(blue and green lines), and spin wave theory^ (thin black 
line). For small J2/J1, the spins configurations are stag- 
gered, while for large J2/J1, the ground state breaks lattice 
symmetry to develop collinear order as shown in the figure. 
SP(N) (bold red line) tends to overstabilize the long range 
ordered phases, most dramatically on the one sublattice side, 
where the critical spin is independent of the strength J2 of 
the frustrating diagonal bonds^. Symplectic- TV restores the 
frustration-induced fluctuations by treating both ferromag- 
netic and antiferromagnetic bonds, on equal footing, which 
corrects this overstabilization. The physical spin, S = 1/2 is 
indicated by a horizontal dashed line. 



At high temperatures and large J2/J1, the first bonds 
to develop are the diagonal bonds, A<j. From the mean 
field equations (JHOJ) and 



d 2 F 

ax 3 " 

d 2 F 
d\dh a 
d 2 F 



d\dh a 
d 2 F 
Ob? 
d 2 F 



dXdA a dh a OA b 



d 2 F 
dXdA a 

d 2 F 
dh a dA b 

d 2 F 



(80) 



where this is a schematic of the seven by seven Hessian 



with respect to A, h x , h y , hd, A 2 



Ay, Arf 



When det% 



changes sign, the decoupled solution is changing from a 
free energy minimum to a maximum, indicating the pres- 
ence of a second order phase transition. By examining 
the unstable eigenvectors, we know which bond fields are 
turning on, without having to solve the seven mean field 
equations. 

All of the matrix elements have similar forms, for ex- 
ample 



d 2 F 
dhl 



d 2 k 1 d 2 u k 



n k (n k + 1) / du> k \ _J_ 
T \dhj J x 
(81) 

Since A and A^ are independent of J\ , we can fix J2 = 1 , 
5 = 1/2 and easily evaluate \ f° r a U Ji & t a given T, 
since the integrals are all independent of J\ . det x = 
can then be solved for J\ c and the phase transition, 
T c mapped out parametrically, as shown in Fig [7] The 
unstable eigenvector is 



—h a 
A„ 



(82) 



showing that the system does develop long range Ising 
order. This method finds all possible second order phase 
transitions, however, it is blind to first order phase tran- 
sitions. As we see in the figure, there is a temperature 
dependent first order transition between the short range 
Neel and decoupled orders which cuts off the second order 
line (see Appendix lAl for derivation). 



1 

T 2 



(rife 



2 LUk 



. 1.2{2c x s v ) 



(77) 
(78) 



we can solve for A and A^ as functions of temperature 
and spin. They are both independent of J\ . A^ turns on 
at a temperature 



T,, 



21og(l + l/5) 



(79) 



Now that we have a full description of the decoupled 
phase, we can look for the next bond fields to turn on as 
we lower the temperature. For simplicity, we assume that 
the spin is large enough that the ground state is the long 
range ordered collinear state, so we know that h x and 
A y must turn on at some point. However, we can look 
for all possible bonds at once by examining the unstable 



1. Analytical form of Trvb 

Now we derive the analytical form for the Ising tran- 
sition temperature, in the limit of large J<ij J\. At tem- 
peratures far below the development of decoupled order 
Xd, but above the Ising transition, T c , the gap in the 
spectrum at (0, 7r/2) 



A sap = VA 2 - (4A d ) 



(83) 



is much smaller than T, and, assuming large 5, we can 
apply spin wave theory to this problem, which implies 
A w 4Ad s» c sw = 4J 2 S. x can be restricted to the 



two relevant parameters h x 
quantities A\, A2 and B 




and A„, and we define the 



(84) 
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FIG. 7: (Color online)Finite temperature phase diagram for 
S — 1/2. Trffequation l79|) and Tn (equation I A2 |) are the tran- 
sitions into short range two sublattice and Neel antiferromag- 
netic order, respectively. The Ising transition, T c is shown for 
both symplectic-TV(blue) and S'P(iV) (red). The Ising order 
is long range, even though the underlying antiferromagnetic 
order is not. The dashed(green) line indicates a first order 
transition from Ising order to short range antiferromagnetic 
order. Just as we saw by examining S c , SP(N) overstabilizcs 
the Ising order. Insets show the appropriate valence bond 
order. 



In the limit A gap — > 0, we find that A\ = A 2 = A = 
B to all divergent orders. This is because our singlet 
bond fields are decoupled from the S = 1 spin waves 
becoming gapless. To find T c , we need to consider the 
short wavelength behavior which makes A — B nonzero. 



detx = (A + B){A -B)- l/J 2 = 



(85) 



A + B is of the order T/A 2 ap , but the divergences cancel 
from A — B and we can calculate this integral to zeroth 
order in A gap : 

1 (d 2 F d 2 F\ d 2 F 
A-B = - ' 



2\dh 2 x <9A 2 J dh x dA y 

~ r rPk ens 

2A 



2 1 d 2 k cos 2 k x cos 4 k y fn k (n k + 1) n k + 5 



(2tt) 2 u 2 V T 

1 f d 2 k cos 2 k x cos 4 k y 
~ ~3TJ (2tt) 2 1 - cos 2 k x sin 2 k y = 

where 7 = .039. Altogether (|85|) gives us 
87 1 

7\2 — T 2 

^gap J l 



7T7 

~T~' 



(86) 



(87) 



We can expand the constraint eauation ([77|) to find the 
gap 



i„„„ / — 8nJ 2 S x 

= exp 



T 



(88) 



which, combined with ([86)1 leads us to the Ising transition 
temperature 

4irJ 2 S 2 



log 



2J 2 S 



while Chandra, 
semiclassicallj^i 



Coleman and 
4irJ 2 S 2 



(89) 



Larkin found 





2J 2 


log 


Jl V 2 lT _ 



(90) 



with 7t = .318. Note that the form of the two temper- 
atures is identical, with only numerical differences inside 
the logarithm, which are negligible for small spin. This 
temperature dependence has been confirmed by classi- 
cal Monte Carlo^i, and quantum numerical studies have 
show that finite S systems also share the temperature 
dependence^. 

The same calculation is much simpler in SP(N) where 
X is a one dimensional matrix, d 2 F/dA 2 ~ — T / A 2 
1/Ji, giving the defining condition j^^PiN) 
1sp(n)I J\- Again inserting the gap (|87p . we find an im- 
plicit equation for j , c SP ( Ar ^ 

rpSP(N) _ 



leTrJaS' 2 



2 
gap 

2 

gap 



log 



16 JgS 2 

T t^S-P(N) 



(91) 



The extra y c SP ( Ar ) j n the logarithm acts to increase the 
Ising temperature, as was also seen in our numerical cal- 
culation(Fig[7|. 

This phase transition has several possible experimental 
realizations. First, there is a direct realization of the two 
dimensional J\ — J 2 lattice in Li2VOSi04, where a transi- 
tion to long range collinear order is immediately preceded 
by a lattice distortion from tetragonal to orthorhombic 
symmetry^. 

In the iron arsenides, a Q = (0, 7r) spin density 
wave order develops either coincident with a tetrago- 
nal to orthorhombic structural transition, or slightly 
belo w 46 ' 47 i 48 i 49 . First principles calculations suggest that 
the system can be described by the J\ — J 2 model with 
J1/J2 ~ although whether the magnetism is itin- 

erant or local moment is still controversial. 

Finally, the spin dimer system, BaCuSi2 6 54 contains 
elements of J± — J 2 physics despite being three dimen- 
sional. The alternating layers of dimers are ordered anti- 
ferromagnetically, but decoupled, like the J 2 sublattices, 
while the interlayer couplings are frustrated like 3\ . The 
compound can be thought of as a multilayer J\ — J 2 
model, where the transition to three dimensionality is 
an Ising transition^. 



V. DISCUSSION AND CONCLUSIONS 

We have identified the time reversal of spin as a 
symplectic symmetry and examined the consequences of 
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maintaining this symmetry in the large N limit. In 
order to write a theory of symplectic spins, all inter- 
actions of the unphysical antisymplectic spins must be 
excluded, leading to a unique large N limit which we 
call symplectic- N. In this paper, we have examined the 
bosonic symplectic- N Heisenberg model. The practical 
consquences are to introduce two mean field parameters, 

^ = <^£* & ]<AU> 02) 

where h{j measures the ferromagnetic correlations along 
a bond {ij} and Ay the antiferromagetic correlations, 
and to identify the mean field theory introduced by Cec- 
catto et al2^ for SU(2) as the unique large N limit. Pre- 
vious large N methods had either ferromagnetism or an- 
tiferromagnetism, and the presence of both means that 
symplectic- N can treat both ferromagnetic and antifer- 
romagnetic states. In frustrated antiferromagnets, this 
is especially important because the frustration manifests 
itself through the presence of ferromagnetic correlations 
on antiferromagnetic bonds; in these cases we call h the 
frustration field. Correctly accounting for the price of 
these frustrated bonds is essential in systems with many 
competing states close in energy. 




b. 



FIG. 8: (Color online) a. The antiferromagnetic triangular 
plaquette has three possible bond orderings: (left) just a 
single A connecting two of the spins and leaving the third 
completely free; (middle) A's(blue) and /i's(red, dashed) seg- 
regated; and (right) the uniform state, which is the ground 
state of symplectic- JV.b. The tetrahedral plaquette. When all 
sites are assumed to be equivalent, there are three different 
types of bonds - as shown in black, gray and thin black lines. 

Frustrated bonds will occur whenever there are trian- 
gles containing two or more antiferromagnetic bonds (see 
Fig. OJtop)). In this paper, we studied collinear magnets, 
where the bonds on the triangle are either exclusively 
ferromagnetic or antiferromagnetic, e.g. - h and A do 



not coexist. In other lattices, like the triangular lattice, 
noncollinear states are expected. Certainly the classical 
symplectic- N limit will contain coexisting bonds, as we 
know hij = S Jij cos fy-, and Ay = S'Jy Sin^-, where 
4>ij 7^ or 7r for noncollinear ground states. Whether this 
coexistence persists in the quantum limit is still an open 
question. In a first attempt, we have examined the tri- 
angular plaquette and found the ground state to be the 
uniform, coextant state. However, the lattice case will 
likely be different; the plaquette version of the J\ — J 2 
model also has a uniform ground state, not the broken 
symmetry state found in the lattice. In the tetrahedral 
plaquette, as in SP(N), there is a continuously degener- 
ate ground state manifold 5 . In SP(N), the degeneracy 
is lifted in the lattice, however, the ground state found in 
the SP(N) semiclassical limit is inconsistent with linear 
spin wave theory 5 ^ Given the many competing states, 
it is an interesting open question whether the frustrating 
fields will bring the lattice ground state into agreement 
with spin wave theory. More generally, we would like to 
know if including the price of frustration substantially 
changes the ground states or response for other highly 
frustrated lattices. 

Now we turn to corrections beyond mean field the- 
ory, the 1/N corrections. These will not affect the phase 
boundaries, but can change the nature of the short range 
phases. Sachdev and Read have shown that the 1/N cor- 
rections for SU(N) spins manifest as a gauge field cou- 
pling to the Schwinger bosons 5 -^. In SU(N), this is a U(l) 
gauge field, which in two dimensions contains instantons 
that generate nontrivial Berry phases which enforce the 
discrete nature of valence bonds. For each spin 5, each 
site participates in exactly 2S valence bonds. The ground 
state alternates periodically between spin-Peierls and va- 
lence bond solid phases as 2S'(modz), where z is the co- 
ordination number of the lattice. Sachdev and Read later 
showed that this treatment can be extended to collinear 
states in SP(N), while noncollinear states do not gener- 
ally have instantons 45 . To examine the effects of h, we 
consider the U(l) gauge symmetry in symplectic- N. 

In the large N limit of the J% — J2 model, the valence 
bond fields A develop between "even" and "odd" sites. 
This breaks the local U(l) symmetry associated with bo- 
son conservation at each site down to a global compact 
U(l) symmetry, under which bi — » e bi on the even sub- 
lattice and bi — > e bi on the odd sublattice (correspond- 
ing to the conservation of J2ieeven n i~J2i£odd Tnc 
instanton tunneling configurations considered by Read 
and Sachdev are space-time monopoles in the electric 
field associated with this U(l) field. In fact, the frus- 
tration fields h link sites on the same sublattice, so that 
h is invariant under the global U(l) symmetry, so it does 
not pick up any phase factor when the instanton forms, 
and it does not modify the the phase factors associated 
with instanton formation. In this way, the frustration 
fields do not affect the formation of valence bond solids 
in collinear states. The effect of the frustration fields on 
noncollinear states is however, still an open question. 
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Another way to move beyond the large N limit is to ex- 
amine the variational wavefunctions which are the ground 
state of the large N limit. The wavefunction of a pure 
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10) 



(93) 



APPENDIX A: J L J 2 FIRST ORDER 
TRANSITIONS 



where Ps projects out the unphysical subspace where 
rib 7^ NS, as given in equation (|45|) . When we include 
the effects of the frustrating fields, 



|*) = Psexp 



<-M> [ | 1°)- 

(94) 

The enp a v ^Ij^j creates effective valence bonds 

of all lengths across the system, causing the spins to 
fluctuate coherently - in the case of the Ising transition, 
these coherent fluctuations break lattice symmetry with- 
out long range magnetic order. 

This paper has addressed the bosonic representation of 
interacting symplectic spins, but the principles of sym- 
plectic closure can equally well be applied to fermionic 
models, either in Heisenberg physics, where Ran and Wen 
have used an identical decoupling 2 ^, or Kondo physics, as 
we have done in collaboration with Maxim Dzero^i. In 
the fermionic spin representation, requiring spins that re- 
verse under time reversal also insures that the spins are 
neutral under particle-hole transformations, which leads 
to a local SU(2) gauge symmetry. In bosonic models, this 
gauge symmetry reduces to the U(l) symmetry discussed 
earlier because db\ ! b^_ a = on site due to symmetriza- 
tion. In parallel with our current treatment of both fer- 
romagnetism and antiferromagnetism in the Heisenberg 
model, we are able to treat both the Kondo effect and su- 
perconductivity within the two channel Kondo model 2 ^. 

The next step is to introduce charge fluctuations while 
maintaining the symplectic spin closure. One possibil- 
ity is to introduce the symplectic- N Hubbard operators, 
which can used to construct the t — J and Anderson mod- 
els. These ensure that a hole hopping onto a site and off 
again will generate a symplectic spin flip. In turn, the 
symplectic closure guarantees that the local SU (2) gauge 
symmetry survives to all orders in N, even at finite dop- 
ing, justifying the SU (2) slave boson theory of Wen and 
Lee in a unique large N limit^. The application of this 
approach as a large N framework for the RVB theory 
of superconductivit y 61 ' 62 is a matter of great interest for 
future research. 
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In principle, calculating first order transitions is simple 
- one calculates the parameters A, h a and A a for each of 
the phases from the mean field equations, plugs them into 
the free energy, or ground state energy at zero tempera- 
ture and compares the energies. In practice, it is difficult 
to solve the mean field equations in complicated phases. 
Second order transitions are much easier because some- 
thing is going to zero. For the zero temperature phase 
diagram of the J\ — J2 model, Fig. [7J we know that the 
transition between Neel and collinear long range order is 
first order because the second order lines(between short 
and long range order of the same type) indicate that the 
phases overlap for S > .4. We have calculated the loca- 
tion of the first order line by comparing the energy of the 
long range ordered states, shown in Fig. [H 




FIG. 9: (Color online)Thin lines indicate the second order 
transitions from short range to long range Neel and collinear 
orders. For points within the region of collinear long range 
order, the ground state energies of the two possible orders 
were compared. Where the collinear order is lower, a circular, 
blue dot is placed on the phase diagram; when Neel order is 
lower, the dot is green and square. The black dot indicates the 
classical, second order phase transition, which was calculated 
analytically. The physical spin, S = 1/2 is indicated by the 
dashed line. 
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As the spins become more classical, the mean field 
parameters become more difficult to calculate, but the 
S — > oo point can be calculated analytically using the en- 
ergies from the previous section ([56]) . For the Neel state, 
4>ij = 7r for nearest neighbor bonds, and for diagonal 
bonds, while for the collinear state <f>ij = n on y and 
diagonal bonds and on x bonds. 



E N = -4Ji + 4 J 2 
E c = -AJ 2 



(Al) 



Thus, the classical transition is second order at J\ 
'2.1 1 . just as found for classical SU(2) spins. The same 
calculation can be repeated with Eg P ( N } , with the same 
result. 

At finite temperatures, there is a first order transition 
between the two short range orders. We already have one 
end of the first order line - the zero temperature point, 



and we can calculate the other end, which is a second 
order point where both antiferromagnetic and decoupled 
short range orders give way to a completely disordered 
high temperature phase. We already know the decoupled 
temperature as a function of J9<|79p. and the antiferro- 
magnetic temperature can be similarly found from the 
mean field equations for A and A in the limit of A — > 0, 



T, 



■Mg+1/2) 
2 log (1 + l/S) ■ 



(A2) 



The two temperatures have identical form - the only dif- 
ference being that where T/v has Ji, Td has J%. Thus the 
first order line ends in a second order point at J\ = J%, 
as show in Fig[7J The intermediate line has been extrap- 
olated, but not calculated. 
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